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Abstract 

We have investigated the dynamics of entanglement between two spin- 1/2 qubits that are sub- 
ject to independent kick and Gaussian pulse type external magnetic fields analytically as well as 
numerically. Dyson time ordering effect on the dynamics is found to be important for the sequence 
of kicks. We show that "almost-steady" high entanglement can be created between two initially 
unentangled qubits by using carefully designed kick or pulse sequences. 
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I. INTRODUCTION 



Control and manipulation of entanglement which is a quantifiable resource for quantum 
information tasks such as quantum computing [1J, communication [2] and cryptography [5] 
have been studied along many directions in the last decade [1HT3]. Among these studies, sys- 
tems that are modelled as 1-D Heisenberg chain with qubits of spin- 1/2 particles as the main 
unit are one of the prototypical examples [H] . For such systems the control of entanglement 
between the qubits can be manipulated with the help of various type external magnetic 
fields P-TT3]. In particular, Heule et al. investigated the feasibility of local operator control 
in arrays of interacting qubits modelled as isotropic Heisenberg spin chains [I] . Along similar 
lines, Caneva et al. explored optimal quantum control by appropriate pulses to affect the re- 
quired transformations by numerical Krotov algorithm [5]. Wu et al. showed that one qubit 
gates can be constructed with global magnetic fields and controllable Heisenberg exchange 
interactions [6]. Levy demonstrated a scheme that uses pairs of spin- 1/2 particles to form 
logic qubits and Heisenberg exchange only to produce all gate operations |7j. Malinovsky 
and Sola studied phase control of entanglement in two qubit systems and showed that by 
changing the relative phase of control pulses, one can control entanglement at will j8] . Sadiek 
et al. studied the control and manipulation of entanglement evolution for a two qubit system 
coupled through XYZ Heisenberg interaction influenced by a time- varying external field [9] . 
Wang et al. demonstrated that near perfect entanglement can be obtained by applying a 
magnetic field on a single spin of an isotropic Heisenberg chain of length N [10J. Abliz 
et al. studied the entanglement dynamics for a two-qubit Heisenberg XXZ model effected 
by population relaxation in the presence of various types of magnetic fields and showed 
that it is possible to produce, control and modulate high entanglement with the help of 
time-dependent external fields despite the existence of dissipation [TT] . 

For a general time-dependent external field, time-ordering effects might be important 
and the dynamics cannot be found analytically. Most of the aforementioned studies employ 
numerical methods to investigate the entanglement control [U EJ [&I-FT0] . Although the nu- 
merical methods are fast and reliable, analytic solutions provide a more clear picture of the 
physics behind the dynamics. For a single qubit, the time evolution of populations and co- 
herence under the influence of external field in the form of Gaussian pulse or a delta function 
kick was investigated by Kaplan et al. and Shakov et al. [I5J HE] • The fast pulse or kick is 
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defined based on the relation between the energy splitting of the qubit AE and the duration 
of the pulse r; if AEt << 1 then the pulse is called a kick [To^ ITS] . Many experimental 
and theoretical studies have been carried out to investigate the dynamics of two level quan- 
tum systems under the influence of a single or a series of kicks for quantum gates [J71 [18] , 
NMR [19], excitation of electronic states in molecules [20] . chemical reactions [21] and quan- 
tum computing [22]. However, there is no study of control of entanglement between two 
qubits by using fast pulses, to the best of our knowledge. 

Fast pulses provide an efficient way of full population transfer in a qubit [I5J [TB] and 
because of that are expected to be an important way of controlling the entanglement. From 
this point, in the present study, we consider two qubits with Heisenberg XXX- type interac- 
tion. Each qubit is under the influence of a local time-dependent magnetic field that acts 
in the ^-direction. We consider one, two, three and four kicks as well as Gaussian pulse 
sequences and their effect on the dynamics of entanglement between the qubits. We show 
that entanglement can be controlled by a careful design of the sequence of kicks. 

The organization of this paper is as follows: In Sec. [XTJ we introduce the model and basic 



formalism necessary to solve time evolution exactly. In Sec. Ill, we discuss time ordering 



effect on the dynamics. In Sec. |IV[ Wootters concurrence as an entanglement measure is 
briefly introduced. In Sec. |Vj the analytic entanglement dynamics of kicked qubits in the 
presence of time ordering is discussed by choosing single and multiple (up to four) kicks. In 



Sec. VI, the effect of finite pulse width on the entanglement dynamics is studied numerically 



by choosing a Gaussian pulse or pulse sequence as an external field. We conclude as a 



summary of the important results in Sec. VII 



II. THE MODEL AND BASIC FORMULATION 

In this paper, we consider two Heisenberg XXX coupled qubits in a time-dependent 
external magnetic field acting in the z-direction. The typical time-dependent Hamiltonian 
for this system may be expressed as [9] (we set h = 1): 

H(t) = H + H mt (t), (1) 
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where 



Hi n t(f) 



j E 

i=x,y,z 
2 



(2) 



i=l 



where ex*' 2 (i = x, y, 2) are the usual Pauli spin matrices, J is the qubit-qubit interaction 
strength and B\(t) and B 2 z (t) are the time-dependent magnetic fields acting on qubit 1 and 
2, respectively. It should be noted that the qubit-qubit interaction term in Eq. is given 
by Hq which is constant in time and the time-dependent part of the total Hamiltonian is 
called Hint(t) which describes the qubit-magnetic field interaction and assumed to be a single 
real function of t. 

The most general form of an initial pure state of the two-qubit system is |^(0)) = 



t=i 



dt 



+a 2 (0)|10)- 


f a 3 (0) |01)+a 4 (0) |00), where a^O) (i = 


a*(o)| 2 


= 1, 


then the probability amplitudes evolve 


; to Schrodinger equation as: 


ai{t) 




" J - B T (t) 
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a 3 (t) 
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(I4 (t) 
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a x {t) 




a 2 (t) 




a 3 (t) 




04(t) 



(3) 



where AB(t) = B 2 z (t) - B\{t) and B T (t) = B\(t) + B 2 z (t). The formal solution of Eq. g 
may be written in terms of the time evolution matrix U(t) as 



(4) 



The evolution operator for the general time-dependent Hamiltonian of two qubits is not easy 
to obtain analytically; a number of systematic procedures are obtained in Refs [23] and [21] 
based on dynamical groups of the system when time-ordering is not important. Here the 



ai(t) 




'cn(O)" 


a 2 (t) 


= U(t) 


a 2 (0) 


a 3 (t) 




a 3 (0) 






_a 4 (0) _ 
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time evolution operator U (t) may be expressed as: 

U(t) = fe-'fo 6 ^'* = fe~ i ft( Ao+6int W) d * 



n=0 



\n ft 



H{t n )dt n ... / H{t 2 )dt 2 / #(*i)dti. 



(5) 



The only non-trivial time dependence in U(t) arises from time-dependent H(t) and time 
ordering T. The Dyson time ordering operator T specifies that H(ti)H(tj) is properly 
ordered [EH M> 125] : 



TH(ti)H(tj) = H(ti)H(tj) + 0{tj - U) H(tj), H(U 



(6) 



where 9(tj — tj) is the Heaviside step function whose value is zero if (tj — ti) is negative 
and one if (tj — ti) is positive. It should be noted that time ordering imposes a connection 
between the effects of H(tj) and H(tj) and gives rise to observable, non-local, time ordering 



effects when 



Hit^.Hiti) ^0 psiEZj. 



III. TIME ORDERING 

If one takes T = 1 in Eq. ^ to obtain U(t), the time evolution is said to contain no 
time ordering. So the difference between the result obtained by an exact treatment of T 
in Eq. ^ and T — > 1 is called the effect of time ordering on the dynamics [151 US]- One 
should note that removing time ordering as T — > 1 corresponds to the zeroth order term in 
an eikonal-like, Magnus expansion in the commutator terms [28J. 



A. Limit of no time ordering 



Replacing T with 1 in Eq. <ph, in the Schrodinger picture we have, 



n=0 



-I) 



III 



H(t')dt' 



E 

n=0 



n ! 



H t+ / H int (t')dt' 



n oo 



E 

n=0 



nl 



H + Hi nt 1 1 



-iHt _ r>0 



(7) 



where H int t = f Q H int (t')dt' is the averaged interaction field, and H = H + H int and 
[H , Hi nt \ terms are non-zero. By expanding in powers of [H(t"), H(t')), it is straightforward 



to show that to leading order in Hi nt (t) and Ho, the time ordering effect is given by 

1 , , rt 



U-U°~-^ fdt" f dt' 



H(t"),H(t') 



H 0) Hl t / dt'(t-2t')f(t') 



(8) 



where H int (t') = H® nt f(t'). This leading term disappears if the pulse centroid T& = t/2 
and f(f) is symmetric about Tf~. Furthermore, U — U° vanishes identically in the special 



i.e., 



cases of H int (t') = 0,H int (t') = H int [El US]. Also the commutator H(t"),H(t') 
the time ordering effect) vanishes for B\{t) = B^(t) or J = because by using the total 
Hamiltonian ([T]) we have 



[H(t"),H(t')} = 2iJ ((Bl(t') - Bl(t')) - (Bl(f) - Blif))) {aft - a^J), 



(9) 



which vanishes when either the time-dependent magnetic fields on qubits 1 and 2 are equal 
or the qubit-qubit interaction is neglected. 

In general there is no simple analytic form for the exact result U(t), except for special 
cases [9j [151 EE]- F° r the result without time ordering with averaged magnetic fields Bit = 
J* Bl(t')dt' = a and B 2 z t = f* B 2 z (t')dt' = 0, the time evolution matrix U°(t) in Eq. Q can 
be easily calculated as: 



u°(t) 



-i(H t+H int t) 

yiy* o 

y(u + iv) y(—w + iz) 

y(w + iz) y(u — iv) 

yty* 



(10) 



where 



y 
m 

u 

V 

w 

z 



JJt 



3 i(<*+/9) 



cos (r) 



— ^ — sin (r) 



0. 



2Jt 
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sin (r) , 



(11) 



where T = 'AJH 2 + (a — f3) 2 and a and /3 are called the integrated magnetic strengths 
associated with the magnetic fields acting on qubit 1 and 2, respectively. 

Similarly, the time evolution matrix without time ordering in interaction picture can 
be studied [IHJ [IS] . However in Refs. [15] and [16], it was shown that the occupation 
probabilities for a kicked qubit for the dynamics in the limit T — > 1 depend on the chosen 
picture. On the other hand, the exact results (the results including time ordering) are 
independent of the chosen picture, as expected. Thus in the next sections, we will discuss 
the entanglement dynamics of kicked qubits by using exact time ordered results, and we 



will work in Schrodinger representation. The matrix in Eq. (10) will be used to check the 
correctness of our results, because as mentioned before, when B\(t) = B 2 (t) or J = 0, the 
time ordering effect defined as U K (t) — U°(t) vanishes. 



IV. MEASURE OF ENTANGLEMENT 



For a pair of qubits, Wootters concurrence can be used as a measure of entanglement [29]. 
The concurrence function varies from C = for a separable state to C — 1 for a maximally 
entangled state. To calculate the concurrence function one needs to evaluate the matrix 

R = p(t)(a y ®a y )p*(t)(a y ®a y ), (12) 

where p(t) is the density matrix of the system and p*(t) is its complex conjugate. Then the 
concurrence is defined as 

C(p) = max{0, Ai - A 2 - A 3 - A 4 }, (13) 

where A 1; A 2 , A 3 and A 4 are the positive square roots of the eigenvalues of R in descending 
order. 

It should be noted that due to the discrete symmetry (conservation of parity under 
flipping of the aj, i = x,y,z and j = 1, 2, i.e. when a\ — > —of) of the total Hamiltonian ([!]), 
the states |$) = a 2 1 10) + a 3 |01) and 1^) = a± |11) + a 4 1 00) can never get mixed in time 
due to that symmetry [9], as can be seen from Eq. Thus we consider the time evolution 
of the concurrence of these states individually. The concurrence function for a pure state 
|$(t)) = a 2 (t) 1 10) + a 3 (t) |01) with density matrix p(t) = |$(t)) ($(t)| is given by 

C(p) = max{0, 2 \a 2 (t)a 3 (t)\}. (14) 
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Similarly, for the pure state \^{t)) = ai(t) |11) + (t) 1 00) with density matrix p(t) 
\^(t)) (ty(t)\, the concurrence function reads as 



C(p) = max{0, 2 |o 1 (*)a 4 (t)|}, 
where according to Eq. (Ej), the time-dependent coefficients read 

4 

Oi(t) = ^^(t)a,(0), 
where £4,(t) (z,j = 1,2,3,4) are the matrix elements of 



(15) 



(16) 



V. ENTANGLEMENT DYNAMICS OF KICKED QUBITS 

In this part, we will examine the entanglement dynamics of kicked qubits by taking into 
account the time ordering effects for the initially pure separable |^(0)) = |01) and maximally 
entangled |$(0)) = ^=(|10) + |01)) Bell states. We will work in the Schrodinger picture and 
present analytic expressions for the time evolution operator of the two-qubit system for a 
single kick as well as a positive followed by a negative kick and a sequence of two, three and 
four equally distanced kicks. For all kick sequences we consider two integrated magnetic 
strength regimes: a = 2(3 and a = 3/3 and for convenience we shall set J = 1 and = 1. 



From those results we will use Eqs. (14) and (16) to analyze and discuss the time evolution 



of the entanglement between the qubits. 

We have also considered the other Bell states ^(|11) ± 1 00) ) and ^(|10) — |01)) as 
well as separable states of two qubits as |11), |00), |10). Under the influence of H(t) of 
Eq. ([I]), the states of the type 1 00) and |11) remain separable, while the entanglement of 



1 11) ± 1 00) ) does not change with time. These may be checked using Eq. (15) and the 



solution of the expansion coefficients a\(t) and a±{t) in Eq. ^ for the considered initial 
states. The dynamics of ^(|10) - |01>) is same as that of ^(|10) + |01)) and that of 1 01) is 
same as 1 10) that are noted after specifying the propagators for kicked qubits and by using 
Eq. dlih and Eq. (fieb . So, we consider only |tf (0)) = |01) and |$(0)) = ^(|10) + |01>) as 



initial states. 

One point we want to emphasize that before the field is active, the propagator is equal 



to e lHot and given by Eq. (23). Based on Eqs. (14), (16) and (23), the concurrence of the 



initial state |$(0)) = 4§(|10) + 1 01) ) is equal to 1, while the concurrence of |^(0)) = 1 01) 
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is | sin(4J£)| at time t before the kick. Note that the entanglement dynamics for the initial 
state |$(0)) = ^=(|10) + 1 01) ) is unperturbed by the qubit-qubit interaction in the absence 
of external field, while qubit-qubit interaction creates a high degree of entanglement that 
oscillates between and 1 if the qubits are initially prepared in |^(0)) = |01) state. 



A. Single kick 




FIG. 1: Concurrence as a function of dimensionless time, Jt, for an ideal positive kick applied at 
Ti = 5 for the initial pure states |$(0)) = -Tg(|10) + |01» (a) and |^(0)) = |01) (b). The dashed 
lines correspond to a = 2/3 and the solid lines to a = 3/3. 



Here we consider two qubits whose states coupled by an interaction field which can be 
expressed as a sudden "kick" at t — Ti, namely B\(t) = a5(t — Ti),B 2 z {t) = (35(t — Ti). 
For such a kick the integration over time is trivial and the time evolution matrix in Eq. ^ 
becomes [To] 

e _j_ff ( t _ Tl ) ^-i H int (t')dt' e -iH T! 



U K (t) 



in the same form as Eq. (10) with elements 

jjt 



V 



Vi = e 

u = cos (2Jt) cos(a — (3), 

v = cos(2J(t-2Ti))sin(a-/3), 

w = sin(2J(t-2T 1 ))sm(o;- / 5), 

z = — sin (2 Jt) cos(ct — (3), 



(17) 



(18) 



for t > T\. The propagator without time ordering is given in Eq. (10) and as explained 
before when a = (3 or J = 0, the time ordering effect, U K (t) — U°(t), vanishes after the field 
is active. 

By inserting Eq. disj) into Eqs. duh and ([Teh for the initial states |$(0)) = 73d 10) + |01>) 
and |^(0)) = |01), one can obtain the analytic expressions of the concurrence functions after 
the kick (t > T\). For the maximally entangled state, the concurrence is given as: 

COS) = max{0, |cos 2 (A) + e 8a( *" Tl) sin 2 (A)|}, (19) 

while for |^(0)) = |01) the concurrence after t = T\ can be obtained as 

C(p) = 2 max{0, |A|}, (20) 

where A = (cos(2 Jt) cos(A) - i cos(C) sin(A)) (i cos(A) sin(2 Jt) + sin(C) sin(A)), A = a- (3 
and C = 2J(t-2T 1 ). 

The dynamics of concurrence for the initial Bell state and the separable state for the 



single kick which are given by Eqs. (19) and (20) for t > 7\ are displayed in Fig. 1(a) 
and (b), respectively. The effect of the kick on the entanglement is pronounced for both 
initial states; the concurrence of the Bell state starts oscillating with an amplitude that 
depends on the ratio of the integrated magnetic strength of the external fields on qubit 1 
and 2 (i.e., a and (3, respectively). The effect of the kick on the system initially in separable 
state, 1 01) , is similar with the exception that the concurrence variation amplitudes get lower 
after the kick. One should also note that C(p) of the initial Bell state is independent of 
J before the kick, while the frequency of its time-dependence after the kick is proportional 
to the qubit-qubit interaction strength J, as can be seen from Fig. 1(a) and the analytic 



expression Eq. (18). 



B. A positive followed by a negative kick 

The propagator for a sequence of either identical or non-identical kicks can be easily 



obtained by multiplication of several matrices of the form of Eq. ( 17) [151 US]- F° r example, 
one may consider a sequence of two kicks of opposite sign at times t — T\ and t = T 2 , 
namely, B]{t) = a (8(t - T x ) - 5(t - T 2 )), B 2 z (t) = (3 (8(t - T x ) - 5(t - T 2 )). Following the 
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procedure given in Eq. (17), one obtains the time evolution matrix for t > T 2 as [T5] : 



U K (t) = e~ ifloi - t ~ T2) e~* ^2-< Hint{t')dt' e _ i ^ (T 2 -T 1 ) e -i H int {t')dt' e -iH 



Ti 



(21) 



where the elements of Eq. (21) is the same form as Eq. (10) with parameters 



y = e u \ 

yi = i, 

u = cos (2Jt) cos(A) 2 + cos (2 J(t - 2T S )) sin(A) 2 , 

v = (cos(Ci) — cos(£ 2 )) sin(A) cos(A), 

w = (sin(Ci) — sin(£ 2 )) sin(A) cos(A), 

z = — sin (2Jt) cos(A) 2 - sin (2 J(t - 2T S )) sin(A) 2 , 



(22) 



where Q = 2J(t — 2Tj), A = (a — j5) and T s = T 2 — T±. For this case, the time evolution 
matrix without time ordering is given by 



U°(t) 



-iH t 



-iJt 



< - 

e a *cos(2Jt) -ie iJt sin(2 Jt) 
-ie iJt sin(2 Jt) e iJt cos(2Jt) 















-iJt 



(23) 



because for a positive kick followed by a negative kick the averaged interaction Hamiltonian, 
Hintt = 0. For the cases J = 0, or T s = 0, or a = (3, the time ordering effect defined as 
U K (t) — U°(t) goes to zero, as expected. 

The entanglement dynamics under the positive-negative kick sequence at times t > T 2 is 



obtained by using the expression Eq. (22) in Eqs. (14) and (16) and are displayed in Fig. 2(a) 



and (b) for the initial Bell and separable states, respectively. The effect of the negative 
kick at T 2 is found to be opposite for the |$(0)) = ^(|10) + |01)) and |^(0)) = |01) initial 
states; for the 1 01) state, the dependence of concurrence on the integrated magnetic strength 
vanishes while for the initial Bell state, amplitude of concurrence oscillations change with a 
and p. One peculiar result from Fig. 2(a) is the observation that the negative kick has no 
influence on the dynamics of concurrence for a = 3/3 magnetic fields (solid line in Fig. 2(a)). 
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FIG. 2: Concurrence as a function of Jt for an ideal positive kick applied at T\ = 5 followed 
by an ideal negative kick at T 2 = 10 for the initial pure states |$(0)) = ^(|10) + 1 01) ) (a) and 
1^(0)) = 1 01) (b). The dashed lines correspond to a = 2/3 and the solid lines to a = 3/3. 



The positive-negative kick sequence also demonstrates the strong effect of time ordering on 
the entanglement dynamics. As mentioned before, based on the propagator for T — > 1 case 



given by Eq. (23), the concurrence for Bell state is always equal to 1, while the concurrence 



for initially separable state has oscillations between and 1 for the times t > T 2 . On the 
other hand, as can be seen from Figs. 2(a) and 2(b), the time-ordered propagator leads to 
different results in the concurrence for both initial states. 



C. Two positive kicks 



To show the difference between positive and negative kicks applied after the first positive 
kick on the entanglement dynamics of two qubits, one may consider a sequence of two positive 
kicks applied at times t = T x and t = T 2 , namely, B\{t) = a (6{t — 7\) + 5(t — T 2 j) , -Bf (i) = 



(3 (5(t — Ti) + 5(t — T 2 )). Following the procedure given in Eq. (17), one obtains the time 



evolution matrix Eq. (21) for t > T 2 in the form as Eq. (10) with parameters 



V 

yi 

u 

V 

w 

z 



JJt 



3 2i(a+/3) 



cos (2 Jt) cos(A) 2 - cos (2 J(t - 2T S )) sin(A) 2 , 

(cos(Ci) + cos(C2)) sin(A) cos(A), 

(sin(Ci) + sin(C 2 )) sin(A) cos(A), 

- sin (2 Jt) cos(A) 2 + sin (2J(t - 2T S )) sin(A) 2 , 



(24) 
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where Q = 2J(t — 2Tj), A = (a — (3) and T s = T 2 — T\. Here the propagator without time 



ordering can be calculated by replacing B\t — > 2a and B\t — > 2/3 in Eq. (10) and note for 



a = (3 or J = 0, U K (t) — U°(t) vanishes, as expected according to Eqs. (pi) and (J9J). 




FIG. 3: Concurrence as a function of Jt for a sequence of two ideal positive kicks applied at T\ = 5 
and T 2 = 10 for the initial pure states |$(0)) = -^(|10) + |01)) (a) and |^(0)) = |01) (b). The 
dashed lines correspond to a = 2/3 and the solid lines to a = 3/3. 

The effect of two consecutive positive kicks on the dynamics of concurrence for two qubits 
is displayed in Fig. 3(a) and 3(b) for the initial Bell state, ^=(|10)+ 101)) and separable state, 
|01), respectively. Comparing the analytic expressions of the time-evolution operators for 



positive-negative and positive-positive kick sequences of Eq. (22) and Eq. (24), respectively, 



along with the Fig. 2 and Fig. 3, the effect of the sign of the kicks in the sequence is to 
change the amplitude of the concurrence oscillations. The oscillation amplitude of C(p) for 
the Bell state as well as separable state increases for the positive-positive sequence compared 
to that of positive-negative sequence of kicks. Also, a/ (3 dependence of the amplitude is 
different as can be seen from a comparison of Fig. 2 and 3. 

D. Three and four positive kicks 



One may consider a sequence of n-positive kicks applied at times t = T\,t = T 2 , t = T n , 

n n 

namely B\(t) = aS(t — Tj), B 2 z (t) = {35(t — Tj). For example, following the procedure 

i=l i=l 

given in Eq. (fT7j), one obtains the time evolution matrix for three positive kicks at times 
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FIG. 4: Concurrence as a function of dimensionless time, Jt, for 4-successive ideal positive kicks 
for the initial pure states |$(0)) = 771 (| 10) + 1 01) ) (a) and |^(0)) = |01) (b). Here the dashed lines 
correspond to a = 2/3 and the solid lines to a = 3/3 and we take T\ = 5,T2 = 10, T3 = 15 and 
T 4 = 20. 




FIG. 5: (Colour online) The contour plot of concurrence versus Jt and the ratio, a//3, for the initial 
pure states |$(0)) = (| 10) + 1 01) ) (a) and |\£(0)) = |01) (b). Here the contour plots include four 
ideal positive kicks applied at T\ = 5, T2 = 10, T3 = 15 and = 20. (There are ten equidistant 
contours of concurrence in the plots between (black) and 1 (white).) 



t > Tq as: 



U K (t) = e ~ iflo{t ~ Tz) e~ l Urnt{t')dt' e -iH { Tz -T 2 ) e -i S^-l H int (t')dt' 



e -iH {T 2 -T 1 ) e -if^H int (l/)dt' e -iH T 1 ^ ^5) 



14 



in the same form as Eq. (10) with parameters 



y = e 



Vi = e 



II 



cos (2 Jt) cos(A) 3 - cos ( 2J (t + 2 ( T i ~ T j))) cos(A) sin(A) 2 



i<j 



W 



Y cos(O) sin(A) cos(A) 2 - cos (2 J(t - 2(T X - T 2 + T 3 ))) sin(A)'- 

i=l 
3 

sin(Ci) sin(A) cos(A) 2 - sin (2J(t - 2(T\ - T 2 + T 3 ))) sin(A) 3 

i=i 

3 

- sin (2 Jt) cos(A) 3 + ^ sin (2J(t + 2(T t - Tj))) cos(A) sin(A) 2 . 



i<j 



Similarly, the time evolution matrix for four positive kicks at times t > T4 

K /+\ -iHoit-n) e -i St^-I H mt (t')dt' e -iH (T 4 -T 3 ) e -i Jt^ H int (t')dt' e -iH (T 3 -T 2 ) 



U K (t) = e 

x e 



i/^+e H int {t')dt' iff (T 2 -T!) -if^l' H lnt (t')dt' iH oTl 



with parameters specified in Eq. (10) 



y = e 



Vi = e 



4i(a+/3) 



cos (2 Jt) cos(A) 4 - cos (2J(t + 2(T, - 1}))) cos(A) 2 sin(A) ; 



i<jr' 

+ cos(2J(t + 2T 123 4))sin(A) 4 , 
4 4 



(26) 



(27) 



^cos(0)sin(A)cos(A) 3 - ^ cos (2 J(t - 2(T - Tj + T k ))) cos(A) sin(A) ; 



4 



i,j,k=l 
i<j<k 



W 



Y sin(O) sin(A) cos(A) 3 - ^ sin (2 J(t - 2(2} - Tj + T fc ))) cos(A) sin(A) J 



i=l 



i,j,k=l 
i<j<k 



sin (2 Jt) cos(A) 4 + ^ sin (2J(t + 2(T; - Tj))) cos(A) 2 sin(A) 2 



i<3 



- sin(2J(t + 2T 12 3 4 ))sin(A) z 



(28) 
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where Q = 2J{t — 27*), A = (a — (3) and T1234 = (Ti — T 2 + T 3 — T 4 ). Here the propagator 
without time ordering can be calculated by replacing B\t — > na and B 2 z t — y n/3 (here n = 3 



for three positive kicks and n = 4 for four positive kicks) in Eq. (10) and note that for 
a = (3 or J = 0, the time ordering effect defined as U K (t) — U°(t) disappears after the field 
is active. 

The dynamics of C(p) under three and four positive kick sequence is shown in Figs. 4(a) 
and 4(b) for ^(|10) + 1 01 ) ) and |01) initial states, respectively. The most important finding 
from these figures is that almost constant high entanglement can be obtained after the 3 rd 
kick for the ^(|10) + 1 01 } ) state and after the I s * and 4 th kicks for the 1 01 ) state at external 
magnetic field ratio of a/ (5 = 2. 

The integrated magnetic field strength dependence of the concurrence dynamics is shown 
in Fig. 5 where we display the contour plot of C(p) as functions of a//3 and Jt for 
^(|10) + |01>) and |01) initial states. For the ^(|10) + |01>) initial Bell state, the max- 
imally entangled state is found to be unaffected by the external field for a//3 — 1 and 
a/ (3 = 2.5,4.25,5.75,7.25,8.75. The a/ (3 periodicity of the maximum of C(p) increases 
after each kick for the initial Bell state. The dependence of C(p) on a/ f3 for the initial sep- 
arable state 1 01 ) is more complicated compared to the case of initial Bell state. The almost 
periodic structures exist also in Fig. 5(b); their periodicity changes after each kick, but it is 
not easy to obtain an expression for that change. Most importantly, the high entanglement 
regions, which are indicated in white in the contour plots have long life times for each pos- 
itive kicks for the initial Bell state, while for the separable state they are distributed in a 
narrower area compared to the initial Bell state case and have long lifetimes only after 1 st , 
2 nd and 4 th kicks. 



VI. ENTANGLEMENT DYNAMICS OF QUBITS PERTURBED BY A SE- 
QUENCE OF GAUSSIAN PULSES 

Depending on the physical implementation of the qubit, it might be difficult to ob- 
tain an external field that can be considered as a kick. Instead a Gaussian pulse with 
finite width can be applied (for example, half-cycle electromagnetic pulses with width near 
t = lps may be experimentally achievable [301 EI])- Thus in this part, we will discuss the 
entanglement dynamics of two qubits under the influence of Gaussian pulses of the form 
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_ (t-Tfe) 

B l z {t) = -^e ^ (ai j2 = a, p) centered at Tt with width r. The dynamics of entangle- 
ment in the presence of time ordering for the initial pure states |$(0)) = ^(|10) + |01)) 
and ^(O)) = 1 01 ) are evaluated by numerically integrating the corresponding equations in 
Eq. Q and using Eq. p4p. Here we will investigate how the entanglement depends on the 



pulse width r by choosing a single pulse, a positive pulse followed by a negative pulse, and 
multiple positive pulses up to four centered at times T% = 5,T 2 = 10, T 3 = 15 and T 4 = 20. 
For all pulse sequences we will consider two integrated magnetic strength regimes: a = 2/3 
and a = 3/3 and for convenience we shall set J = 1 and (3 = 1. One should note that in the 
limit t — > 0, the results of entanglement dynamics of kicked qubits in the presence of time 
ordering should be the same that are analyzed in the previous section. 

A. Single Pulse 

In Fig. 6, we show the results of a calculation of the concurrence for the initial pure states 
|$(0)) = ^(|10) + |01)) and |$(0)) = |01) when strongly perturbed by a single Gaussian 
pulse centered at t = T\ with width r. According to Schrodinger equation ^ of the system 
considered here, the expansion coefficients a\(t) and a^(t) evolve independently and a 2 (t) 
and a^it) obey a first-order coupled differential equation set, for example for a single pulse, 
as: 

ia 2 {t) = I — J — v ^ e a 2 (t) + 2Ja 3 (t), 

zd 3 (t) = (- J + ^^-e-^^j a 3 (t) + 2Ja 2 (t), (29) 

which are solved numerically by using a fourth order Runge-Kutta algorithm. The most 
important observation from Fig. 6 is the existence of almost constant high concurrence for 
the initially separable state at a = 2/3 integrated magnetic strength and high width Gaussian 
pulse, while the entanglement continues to have high amplitude oscillations for a = 3/3; its 
value for a = 2/3 is almost constant at around 1 for Jr; the dimensionless pulse width 
greater than 0.15. On the contrary for the initial Bell state, the oscillation amplitude of 
C(p) increases with the Jr of the pulse for each magnetic ratio (ct//3 = 2 and a//3 = 3). 
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fcj.jE=lI.05 0-') Jt=u lis 




FIG. 6: Concurrence as a function of Ji for a single Gaussian pulse with width r for the initial 
pure states |$(0)> = -^(|10> + |01>) (a), (c), (e) and (g) and |tf(0)) = |01) (b), (d), (f) and (h). 
The dashed lines correspond to a = 2/3 and the solid lines to a = 3/3. Here we assume four 
dimensionless pulse width as: Jr = 0.05, 0.1, 0.15, 0.2. 

B. Positive- negative and positive-positive pulse sequence 

In Fig. 7 and 8, we show the results of a calculation for the concurrence for the initial 
pure states |$(0)) = 4=(|10) + 1 01) ) and |^(0)) = |01) when strongly perturbed by a single 
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Jt 



FIG. 7: Concurrence as a function of Jt for a positive followed by a negative Gaussian pulses 
having the same width r for the initial pure states |3>(0)) = (| 10) + 1 01) ) (a), (c), (e) and (g) 
and 1^(0)) = 1 01) (b), (d), (f) and (h). The dashed lines correspond to a = 2(3 and the solid lines 
to a = 3(3. Here we assume four dimensionless pulse width as: Jt = 0.05, 0.1, 0.15, 0.2. 

Gaussian pulse centered at t = T\ followed by a negative or positive Gaussian pulse centered 
at t = T 2 with the same width r. For the double pulse sequence a 2 (t) and a 3 (t) obey the 
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coupled equations: 

za 3 (t) = f - J - (e ^±e ^~)J a 2 (t) + 2 Ja 3 (t), 

ia 3 (t) = \^-J+ ^ (e ^±e J a 3 (t) + 2Ja 2 (t), (30) 

where + sign in the ± on the right-hand side is for positive-positive pulse sequence, while 
— sign is for positive-negative pulse sequence. 

The concurrence dynamics for positive-negative and positive-positive pulse sequences 



obtained from the numerical solutions of Eq. (30) are displayed in Figs. 7 and 8, respec- 
tively for the initial Bell and the separable states at different dimensionless pulse width 
values (Jr = 0.05,0.10,0.15,0.20). For the initial Bell state under positive-negative se- 
quence, the most important effect of the pulse width seems to be an increase in the C(p) 
oscillation amplitudes for a = 3/3 at times t > T 2 . The almost constant high entanglement 
can be obtained for the 1 01} initial state for positive- negative pulse sequence as can be seen 
from Fig. 7(h) for a = 3/3 and t > T 2 . One peculiarity of this figure is that high entan- 
glement is obtained for a = 2/3 after the first pulse, while it is obtained for a = 3/3 after 
the second pulse. For the positive-positive Gaussian pulse sequence, the difference from 
positive-negative sequence becomes small as the width of the pulse gets larger as can be 
deduced from a comparison of Figs. 7 and 8. On the other hand, for a small pulse width, 
the difference is significant. For example for Jr = 0.05 and a = 3/3, the initial Bell state 
has nearly constant entanglement around 1 (see Fig. 7(a)) for positive-negative pulse se- 
quence after negative pulse, while it oscillates between 1 and 0.25 for positive-positive pulse 
sequence (see Fig. 8(a)). 



C. A sequence of four positive pulses 

The effect of integrated magnetic strength and the pulse width on the dynamics of con- 
currence for two qubits perturbed by a sequence of four positive Gaussian pulses is dis- 
played in Fig. 9(a)-(h) for initial Bell state, |$(0)) = ^(|10) + |01)), and separable state, 
1^(0)) = 1 01). For this four positive pulse sequence the concurrence may be calculated by 
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FIG. 8: Concurrence as a function of Jt for a sequence of two positive Gaussian pulses having 
the same witdh r for the initial pure states |$(0)) = ^(| 10) + |01)) (a), (c), (e) and (g) and 
1^(0)) = 1 01) (b), (d), (f) and (h). The dashed lines correspond to a = 2/3 and the solid lines to 
a = 3/3. Here we assume four dimensionless pulse width as: Jt = 0.05,0.1,0.15,0.2. 
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using the numerical solutions of the coupled equations: 



ia 2 (t) = 

ia 3 (t) = [-J + 




a 3 (t) + 2Ja 2 (t). 



a 2 (t) + 2Ja 3 (t), 



(31) 



These figures can be compared with Figs. 5(a) and 5(b) to discern the effect of the pulse 
width. As the pulse gets wider, the a//3 dependent oscillatory structures in the figure coa- 
lesce to produce non-periodic structures, especially after third and fourth pulses. The high 
entanglement regions, which are shown in white, still can have long lifetimes, as indicated by 
white straight perpendicular sections in the contour plots of Fig. 9. In the case of ideal kick, 
the maximally entangled state is found to be unaffected by the external field for a/{3 — 1 and 
a/0 = 2.5,4.25,5.75,7.25,8.75 (see Fig. 5(a)). Comparing Fig. 5(a) with Figs. 9(a), (c), (e) 
and (g), the initial Bell state is found to be unperturbed by the highly wider Gaussian pulses 
if and only if a/0 = 1; especially seen obviously for the dimensionless pulse width greater 
than 0.15. One should note that this is one of the conditions in which time ordering effects 
vanishes. 

The last point we want to emphasize that what happens the entanglement dynamics 
between two qubits if the time ordering effect vanishes. As mentioned before the time 
ordering effect vanishes for the special cases either a = or J = 0. From corresponding 
equations it can be noted that for the case a = 0, the concurrence function for the initial 
state |$(0)) = -^(|10) + 1 10)) is equal to 1, while for the initial state |^(0)) = |01) is equal 
to | sin (4Jt)| and unaffected by the sequence of the kicks or Gaussian pulses. For the other 
case, J = 0, under the influence of kick or Gaussian pulse sequences the concurrence for the 
initial Bell state is always equal to 1, while the separable state remains separable at any time. 
These show an important fact that since the time ordering provides a connection between 
interactions at different times, it is responsible from the nonlocal correlations between the 
qubits in time. 



VII. CONCLUSION 



We have investigated the dynamics of entanglement for two qubits that interact with 
each other via Heisenberg XXX-type interaction under a time-dependent external magnetic 
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field. Initial state of the system is considered to be pure Bell or separable states. The main 
aim of the study was to investigate the controllability of the entanglement with a sequence 
of pulse or kick type external fields. 

The effect of time ordering in the dynamics of concurrence is found to be important; 
concurrence calculated when the time ordering is neglected is found to be completely different 
than when it is taken into account. Time-dependent concurrence obtained after one, two, 
three and four kicks at different magnetic field strengths indicate that one can employ 
carefully chosen kick or kick sequences to produce high entanglement between two initially 
non-entangled qubits. 

We have also considered the effect of the pulse width of the external field on the en- 
tanglement dynamics by modelling the external field as a Gaussian pulse or a sequence of 
Gaussian pulses. Increasing the width of the pulse is found to enhance the control of high 
and steady entanglement. 

One should note that the external control field considered in the present study acts on 
both of the qubits at the same time. It might be possible to use pulse sequences acting an 
individual qubits at different times to obtain a better control of entanglement. 
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FIG. 9: (Colour online) The contour plot of concurrence versus Jt and a//3, for a sequence of four 
positive Gaussian pulses of width r for the initial pure states |3>(0)) = ^(|10) + 1 01} ) (a), (c), (e) 
and (g) and |^(0)) = 1 01) (b), (d), (f) and (h). Here we assume four dimensionless pulse width as: 
Jt = 0.05, 0.1, 0.15, 0.2. (There are ten equidistant contours of concurrence in the plots between 
(black) and 1 (white).) 
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